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ABSTRACT 


The system studied is a surface-to-air guided missile whose 
aerodynamic and control characteristics can be represented by a set 
of linear differential equations. The objects of the study are to select 
a set of control system gains which will give specified performance af 
the missile within the roll and yaw subsystems and at the same time 
to minimize the coupling effects whereby motion about one of the axes 


causes dynamic response about the other axis. 


The system is studied analytically using the mathematics of 
control system synthesis and design. Analytic investigations are 


paralleled by system simulation on a REAC analogue computer. 


Using the above methods, gains are selected which give 
specified response in the roll stabilization and yaw control subsystems 
and suitably damped response due to coupled motions between systems. 
The study also investigates the possibilities of compensation to 


attenuate or eliminate coupled response. 
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OBJECT 


The object of this paper is to analyze a surface-to-air missile 
using equations which describe its aerodynamic and control system 
responses; to select values for control system parameters which will 
give specified response; and to investigate the possibility of compen+ 
sation which will eliminate coupled response effects between sub- 


systems. 





CHAPTER 1 


INTRODUC TION 


ee Physical Description of the Missile System 


The missile studied is a surface launched antiaircraft weapon. 
It is a homing missile using proportional navigation. Thus the missile 
borne receiving antenna tracks the target and the guidance system 
continually measures the angular rate of change of the line of sight to 
the target. The guidance computer generates command signals for 
yaw and pitch autopilot systems. The autopilot systems react by 
changing the missile heading at a rate proportional to the rate of change 
of the line of sight. Because of cruciform missile shape, the yaw and 
pitch autopilots have the same form. Polarization of the receiving 
antenna and the insertion of additional autopilot commands in earth 


coordinates require roll stabilization of the missile. 


The general control problems of the autopilot, in response to 


commands generated by the guidance computer, may be stated as: 


We An output acceleration in yaw equal to the command 


acceleration in yaw with as fast a response as possible, 
Ze Same response as above, but in pitch. 


OR Maintenance of roll stabilization (zero roll angle) in the 


presence of induced roll torques. 


1,2 System Equations 


The coordinate system to be used is shown in Figure l.1. 


The axes are centered at c.g. and rotate with the missile airframe. 





Fig, 1.1 Missile Axis Orientation 





Consider a 7, % to be unit vectors along the x,y,z axes re- 


spectively. Define the following: 
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applied force = Fi + BY + F 
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velocity of the missile = er at wk: 


angular velocity of the coordinate system 


pi + qpt+ rR 
mass of missile 
applied moment = M4 a M,} + M 


angular momentum of missile 
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be principle axes. 


motion are: 


Consider the missile to be a rigid body, and the axes chosen to 


m (u - rv + qw) 
m (Vv - pw + ru) 


m (w - qu + pv) 
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Iq io (1, = I) rp 


Lr - (1, - 1) pq 


The force and moment component equations of 
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If we consider the gravity components of force as inputs to the 


system,which are trimmed out in steady flight by command bias signals 


we may write the applied forces and moments as: 


where: 
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Rocket thrust, aligned with x axis 
Dynamic pressure = =p vie 

Mach number 

Reference area, (usually body cross section) 


Nondimensional drag coefficient 


Nondimensional yaw force coefficient 
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Cy = Nondimensional pitch force coefficient (1, 
C) = Nondimensional rolling moment coefficient (1, 
C4 = Nondimensional pitch moment coefficient Cie 
Oo =  Nondimensional rolling moment coefficient (1. 
d = Reference moment arm (usually missile 

body diameter) cle 
B = Yaw angle of attack = arc sin a ce 
Q = Pitch angle of attack = arc sin = (1. 
1 = y = Control surface deflection to produce 

yaw forces and moments (1, 

= = Control surface deflection to produce pitch 

forces and moments ce 
6 = Control surface deflection to produce rolling 

moment le 


From Figure 1.2 we see that y , 1, 


actual control surface deflections ® by 
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Fig, 1.2 Positive Deflection of Control Surfaces 


Now consider the missile to be in the following steady state 


condition: 
1. Constant speed (constant Mach) where 


reese? 2 2 _ f. 2 2 
V=Vo= Nut vot wi = Vu + wy, (1. 38) 


2, Missile trimmed to constant pitch angle of attack, a, 


This means that 


M 
iy 


QSdC_ = 0 (1, 39) 


and EF 
Z, 


Ht 


QSCx = constant (1. 40) 


The above may be expressed by: 


v = p =r =8 =6 Ne = 0 (1. 41) 





and 


i constant (1. 42) 


Dp 


as constant (1. 43) 


Allow motions of v, r and p as small perturbations about the 
steady state values noted above, and consider their products to be 
negligible. We find, on substituting into the equations in which v, r 


and p occur: 


ac. aC, SF 
ee eee tn, Y" a5 
° ; 
= m f[ V5 (cos B) £- wie (sin as) p+ Ve (cos as) Te (1. 44) 
ce a coC oe dC dC 
a pol eal oa. Sel aes - 20 (1, 45) 
M, = Q8d(50 P+ 5g 8 +55 Oo + 94 eee 
aC. oC, oC. . 
M, = QSd (FF B+ pe tt oie) Baa leas pages (1. 46) 


The partial derivatives above are stability derivatives!'+ and 
are obtained empirically for airframe and conditions specified by wind 


tunnel tests. 


Substitute the below assumptions and definitions into equations 
(1.44) - (1. 46) to arrive at the final linearized aerodynamic performance 


equations of the missile, equations (1.64) - (1. 66). 


cos B = 1, since B less than 15° (1. 47) 


@ = constant since V is constant and p is constant 
at steady altitude (1. 48) 


1a; I. and I are constants since propellant burning 


is slow process (1, 49) 

p = @ = roll rate (1. 50) 
3 

ieee e= pitch rate (1; 51) 
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Table A.1, Appendix A, tabulates the above aerodynamic 


coefficients at three pitch angles: 
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Yaw force equation: 


(Ab + By + Mé)= 8 -¢ sina +h cosa (1. 64) 


Yaw moment equation: 


ae I. Sal —«@ Negligible 
Clee at Nic = i= { oo. term (1.65) 


Z 


Roll moment equation: 


Gé+ Fo+HB+ Ly=%$ (1. 66) 


The above development is based on information received from 
Convair (Pomona) and roughly parallels a more rigorous development 
of linearized aircraft performance equations outlined in Reference 1.1, 
Sections JI « I through II - 8. 


Having the aerodynamic performance equations we need now to 
consider the control systems in yaw and roll (pitch being considered 


constant in this analysis). 


The following sensors are used in the missile: 


F 
l, Accelerometers which in effect measure = and = ; 
or with suitable scale factors, (AB + By+M65) and (Aa + By+M6). 
2. Rate gyros which measure $ ; rf : 6 : 


The accelerometers provide information for feedback 
comparison with command signals and the angular rate signals are 


used for damping. 


The sensor outputs used in this analysis are related to the 


sensed quantities as shown below: 


(AB + By+M65)ti(s) = G,(s)(AB + By+M6)(s) (1.67) 
b'(s) = Go(s) d(s) (1. 68) 
§'(s) = Ga(s) b(s) (1. 69) 





where the primed quantities are the sensor outputs and the sensors are 


quadratic of the form 


2 
wh 
G(s) = 5 


2 
s + bw 8 tw 


5 (10) 


Nn 


The control system equations for the autopilots considered in 


this analysis are the yaw control equation 


y,(s) = [ K, (AB + By+ M6) (s)+ (Ky + Kgs) $'(s)} (1.71) 


and the roll control equation 


e Ky e 
6 (s) = - [ a a7 K, + Ke s|¢o'(s) (1. 72) 


These are integrated through rate servos to control surface 








deflections 
G(s) , 
00) eee 5) (1.73) 
and G, . 
6(s)= : 6 (s) (1. 74) 


where G, and G, are of the form indicated in equation (1.70). 


A functional diagram representation of the system is shown in 


Pieure 1,3. 


1.3 Performance Specifications 


Consideration of the equations representing the aerodynamic 
response functions of Figure 1.3 shows that interaction exists between 
(AB + By + M6), ¢ Ju and y and 6 at oe Oe, 


Since the system can be decoupled at as One specifications for 


response at that pitch angle are to be met first. As set forth by 
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Aerodynamic 
Response (AB +B +Mé6) 


Functions 


K4 . 
me +K, + Kes) | 
(AB +By ), 
= / 
(AB+By +Mé6) 





Fig, 1.3 Missile System Functional Diagram 





Convair in stating the problem, the responses to be achieved (in 


uncoupled mode) are: 


ay Time constant of the yaw subsystem to step command 
(AB + BY). to be less than .6 seconds. (Time constant here is defined 


as the time for the quantity to reach 63% of its steady state value. ) 


2. Time constant of the roll subsystem response to a step 
torque (6,7) to be less than .5 second. (The time for ¢ to decay to 
less than 37% of its peak value). 


The only requirement in coupled modes (a, = en 24°) is that 
the system be stable and reasonably damped in response to step 


torque, be , or step command, (AB + Bae : 


1] 





In addition, other considerations have placed the following 


limits on Ke and Ke: 


K, < . 40 (1.75) 
Kz S02 (1, 76) 
All gains (K) oe Ke) are to be considered positive. 


Thus we have the system, and the performance specifications 
to be met. The remainder of this chapter will describe briefly the 


steps to be followed in the system analysis. 


wis 4 Idealized System 


Consider a system in which the servos and sensors have unity 


transfer functions: 


Geiser Ge =l ) 


The missile with "perfect" sensors and servos will hereafter 
be called the idealized system. The idealized system will be con- 
sidered first, and then the effects of non-ideal servos and sensors 
will be discussed. Itis sufficient to say that the idealized system is a 
reasonable approximation to the actual system if it can be shown that 
the servos and sensors introduce into the actual system additional 
modes of low amplitude which are far removed in frequency from the 
natural modes of the idealized system. This matter will be discussed 
in Chapter 6. For the idealized system we may consider the aero+ 


dynamic and control equations in the following form: 


(AB +By+M6)=8-¢ sina + cosa. (1. 78) 
CB +Ey+N6 =U (1.79) 
GS + F¢'+HB+Ly=¢ (1, 80) 
veal K (AB + By+ Mé) + Ky 4 #K, Ul (1, 81) 


12 





§=-—[K,¢+ Ke $ + Keo J (1. 82) 


1,5 Procedure of Investigation 


The investigations of subsequent chapters will 


iF, Discuss analytic methods of selecting gains to provide 


specified performance, 


a Use analogue simulation (guided by the methods discussed 


above) to select a Suitable set of gains. 


3. Discuss the possibility of compensating networks within 
the missile to eliminate or attenuate the effects of coupling on desired 


response. 
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CHAPTER 2 


ANALYTICAL STUDY OF THE COUPLED SYSTEM 


2.1 General 


In Section 1. 4 of Chapter 1 are shown the linearized differential 
equations of the idealized system. These can be put into the form of a 
"system determinant" by gathering all of the variables in each equation 
to one side, the other side being zero. If this be done, and LaPlace 


transform notation be adopted, there results the system determinant: 


% € 
B W od ¥ 6 
(s-A) COS @ -sin a =15 -M 
O O 
-C S 0 -F -N 
KE -H 0 (s-F) a -G (2.1) 
x 
) 0 5 0 l 
S 
-K,A ~Y 0 (s-K,B) -K,M 
h ee eae 
Where = aS + a Tr 4 
ae (K,s + Ko) 


Expansion of this determinant yields the system characteristic 
equation. This is a sixth order polynomial in s in whose coefficients 
appear a considerably involved mixture of the controller gains. 
Appendix B gives this characteristic equation where the aerodynamic 


coefficients for ay 24° have been entered. 
The objectives of the analytical investigation are: 


i To determine how to fix upon a set of gains K) through 


14 





K, which produce a stable system and 


6 
2. To determine a sei of gains which will provide for meeting 
the specified preformance at ae O described in Chapter 1, Section 1. 3 
and 
os To select these gains not only to meet objectives 1. and 2., 
but further to obtain conditions in which the disturbance in roll caused 
by insertion of yaw command and the disturbance in yaw casued by roll 


torque disturbances are acceptably low and of reasonable damp ing. 


Inspection of the characteristic equation reveals that several 
inequalities can be obtained from the expressions for the coefficients 
of powers of s which must be satisfied to prevent occurence of roots 
of the characteristic equation in the right half of the s plane. Given 
a set of gains, the test for presence of sign change in the equation 


can be readily performed. 


It would also be possible to apply the Routh Criterion, at the 
expense of additional labor, for a given set of gains, to obtain further 


information as to the nature of the roots. 


These, and other procedures that might be derived afford little 
in the way of a guide to selection of a new set of gains to test when one 


set has resulted in an unstable condition. 


There is the possibility of forming from (2.1) a set of non- 
interacting transfer functions from which a block diagram of the system 
could be constructed. With manipulation of this diagram and the use of 
root locus methods, the roots for a given set of gains could be found. 
Unfortunately there would be at least three inner loops of such a diagram 
whose roots are gain dependent, and a very considerable amount of 
tedious and repetitious graphical work would be involved, even if each 
new attempt could be based on a set of gains choice of which was 


guided by the last trial. 


It is apparent that it would be desirable to have a more orderly 
approach in which less futile graphic or algebraic work is expended in 
arriving at the conclusion that the last guess was unsatisfactory before 


making anew, more or less intelligent, guess. 


ld 





The object of the following section is to discuss a procedure 
which, while not avoiding trial and error and not avoiding tedium, at 


least puts to good use the experience gained by past results. 


This procedure begins with manipulations of the characteristic 
equations of the uncoupled roll and yaw subsystems. These equations 
are of low order and quite straight forward so that the effect of gain 


changes on the resultant responses are readily seen. 


2.2 An Algebraic Approach to a Stable System Characteristic 
Equation 


In Chapter 3 is derived a block diagram interpretation of the 
system equations which results in a new form of the system determi-+ 
nant where in the coupling effects and the subsystems are clearly 


depicted. 


When (3. 7) of, Chapter 3 is reformed for & = 0 we have 


Yaw Subsystem 


1 -F -F}, | 0 0 
=! l a | 0 0 
= -F, -F) l | 0 Q (2, 2) 
ooo Tm 
QO Q Q | “By il 


Roll Subsystem 


Expansion of (2, 2) will reveal a result which is the product of 
the results obtained by expanding those sections of (2. 2) labelled Yaw 
subsystem and Roll Subsystem as if they were themselves determi+ 
nants. 


When the block diagram Fig. 3.7 of Chapter 3 is considered 
ae, = QO that diagram divides into two new diagrams of the roll and 
yaw subsystems. Using the methods described in Chapter 3, one can 


obtain the subsystem determinants which are the subdeterminants 
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Wathin (2.2). 


Then if we let 


! ag Peer) 1 ein Be ai s + K,) 
2 : 6 5 4 
AF |p 1 | oy, wenn 
ll s (s-F) 
and 
: oe 12 
ay = -F, iL -F, = 
-F, -F, i 
‘3 2 
[q +(-K, B-A-EK,) § + (AK, B-C-EK, + EAK, - BCK, - BK, A) § . 
S(s-K)B)( \-A) 
(CK,B = EAK, + EAK, - BCK,) 
a (s-K BYNs-A) (2, 4) 


Employing the same analog with Cramer's Rule as is used in 
Chapter 3 we can express the outputs due (AB + By). and Orr in this 


ao =0 case as 





, " Fig 
(s) _ Wl l (2 
? a) ara a 
bmp S 4. 
: “Fig | 1 Fig -F ye 
-Fe 1 0 oo Al -F, 
eaeBy(s) “alts “hi VY, BFal® “Fi? (2,6) 
(AB + By(s x A 
y y 


When (2.5), (2.6) are expanded, substitutions made for the defined 
transfer functions, and simplification made, applying a step input 


(AB + By).(s) = ~ and 6,.(s) = + 





BS (2.7) 


o(s) = a: ee mn ea 
(S + (GK, - F)S +K,S + K,) 


“Ky Paes. = i] 


(AB bay )(s)= So a eae a 
Sf. + (-K,B - A - EK,)S + Ste, | 


HS IS [S (S-A)+C] 


+ a. Ol adr co ey ee 
s[s + (-K, B -A- EK.,)S oe | 


In the results of (2.7) and (2.8) it is clear that the modes which 
exist are completely specified by the denominators and that in each 
case these are 3 in number. It is further apparent that the zeros of 


both closed loop functions are due to aerodynamic coefficients. 


ene has derived a large number of forms of ''desired" 
characteristic equations on the assumption that the performance is 
dominated by a pair of complex roots and has prepared graphical aids 
to the rapid determination of the pole-zero configuration of the closed 
loop function which reflects a specified time-domain response toa 
step input. With these it is possible to rapidly obtain a set of gains 
for the separate subsystems to achieve a specified response. However, 
there is no certainty that the assumption of "2nd order dominance" in 
the subsystem will produce a set of gains which will make the coupled 


system stable. 


It is possible, though, to assume a set of gains and rapidly 
determine the residue coefficients of the modes that exist using methods 
described by Thaler. ae By this method a check can be quickly 
obtained upon the ability of a stated set of gains to meet the specifie 
cations at oan QO. A map can then be prepared of the locations of roots 
of the characteristic equations of the subsystems which do satisfy the 


a= O.specifications. 


Having obtained a starting set of values of the gains, these can 
be entered into the coupled system characteristic equation (Appendix B) 


and that equation factored, the roots plotted with the roots of the 
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characteristic equations of the uncoupled subsystems, and the trend 


noted, 


A second trial for the roots can be guided by the concept that 
increased stability of the coupled system will be accompanied by 
increased "sluggishness" in the a. = 0 modes of the subsystems. In 
the actual carrying out of this process it is also noted that transient 
oscillating frequencies of Oe ss Q modes tend higher from the first 
choice of subsystem root locations. By following a path on which the 
complex pair of roots of each subsystem increase in natural frequency 
but remain at constant damping ratio the real root moves toward the 
origin, and convergence to stability of the coupled system characteristic 


equation is apparently more rapid than on other paths. 


Thus there does exist one orderly, albeit tedious, process for 
converging upon a set of gains which produces stability of the coupled 


system and yet allows the subsystems to meet the a, = 0 specifications. 


"Bold spotting" allowed convergence to a stable set of gains in 


three tries using these concepts. 


Somewhat less tedious labor is involved if one does not observe 
the constraint of specifications too rigorously, merely working with the 
subsystem characteristic equations to produce stability in the coupled 
system characteristic equation without too much regard for specifications, 


then trimming to achieve the specifications at a = 0, 


Such a procedure is, then, capable of satisfying objective (1) of 
Section 2.1. It is quite apparent that repetitive factoring of the coupled 
system, characteristic equation as would be necessary with this (and 
most other) methods to reach the other objectives is simply beyond 


reason without automatic computing facilities. 


The conclusion is therefore drawn that this method does set forth 
a concept which is useful in guiding one to a stable set of gains for an 
analogue simulation of the problem, but is otherwise of academic interest, 
and totally devoid of practicality in achieving the other objectives of 


Section 2.1. 


For this reason all work in "trimming up" the system was done by 


i 








analogue simulation as discussed in Chapter 4. 


2.3 Steady State Analysis of the System 


In the preceding section it-was concluded that setting the gains 
in the system to meet all the stated objectives is not a practical task for 
analytical study, that rather it is an undertaking more appropriate to 
the analogue computer. Considerable useful information for a com- 
puter study can be gained by analyzing the features of the variables 


oonsidered as outputs in the steady state. 


The methods described in Chapter 3 for solving the equations 
will be used here. The uncoupled mode at a O will be dealt with 


irst. 


1. Uncoupled Roll Subsystem 


In Section 2.2 it was shown that at a, =0 the roll system charac- 


toristic equation is given by 














I -F 
a s“(s-F) i Gr (1S 3° +K,s + K,) 
a 6 O + 
er a aa (2.3) 
Ss (Ss my F) 
It was also shown that 
5 0 “Fo : 
$ (s)7 Dt (2.5) 
T A ¥ 
I 
F'10 
Since a is analytic in the RHP, the final value theorem is 


| 7 R 
applicable. For bp (s) = ae 


e (= s = F) | 
t> A g > 0 S + (GKe - F)s + GK,s +GK, 
8° (s- F) 


(2.79) 
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Thus the roll system ultimately recovers from a roll step 


disturbance, 


2. Uncoupled Yaw Subsystem 


From (2.4) the characteristic equation of the uncoupled yaw 


subsystem is 


[ s°+(-KB- A- EK,)s° +(AK,B - EK, 
y s(s = K, BXs - A) 


+ EAK, - BCK, - BK,A)s 
A 3 3 ] 


(CK, B - EAK, + EAK, - BCK,)| 
i l y l 2 (2. 4) 
s (s - K,B){s - A) ad 


In (2.8) for a step input (AB+ By) (s) = = 


- PAK, [ Bs - E] 
eS k=) = —_ rr 
Sa lpseesa ics: $< (BC - AE)(K, + K,)] 


- PBK, [ ao -As+C| 
i" . EEE SRA (2, 5) 
| a (BC- AE)(K, + Ko)] 


Then since for a stable choice of gains this function is 


analytic in the RHP the Final Value Theorem applies and 


Lim (AB+ By)(t) = Lim [ s (AB+ By )(s)] 


t>- © s-0Q 
- PK, (BC - AE) K) 
a ae (ee 


Applying the Cramer’s Rule method of solution to the yaw 


subdeterminant in (2, 2) 
-K 
[ Es (s-A)+BC] 


. pe . ale (-F. - F. F.) _ s(s-A)(s-K)B) 
(AB+By)., A, Se (2, 12) 


2] 





or 





° -K, [ Es + (BC - AE)] 
Era ()= nm 
ic MemEse es... (BC - AE)(K, + K,)] 


Using (AB + By) (s) = ~ again, and applying the Final Value 


Theorem 
| “ : -K, J= 
Lim W(t) = Lim s [ w(s)] = K +k. (2213) 
t > 09 s2>0 l Zz 


9 
Thus itis apparent that both (A8+ By)(t) and J(t) come toa 
constant non zero value in steady state and that J(t) increases without 
limif . Their senses and magnitudes are easily calculated for com- 


parison with computer results. 


3. Coupled System 


In performing the steady state analysis at on 0 it is convenient 
to let s ~ 0 in the system determinant (3.7) at the outset since that 
determinant appears in the denominator of the Cramer's Rule solution 
for each variable of interest. In the following statements, the fact that 


this has been done will be signified by [A] - where 








Ba. 4» 
-E 0 -C 0 -N 
Lim A=-~ |= —— l — My-+ifal_. (2.14) 
s>0O 5s Ss 
fF ° - } $ 
me 0 oe 0 


Other notation is as explained in Chapter 3. It is assumed 
throughout that a stable situation exists so that the Final Value Theorem 


is applicable. 
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a. go due to roll disturbances bp (s) apie 


Ss 
gee 
Cees (2,15) 
R 
sl 4,,I ss Bee ee oe 
Lim ¢(t) =Lim s [ ¢(s)] = i = an (2.15a) 
t 00 s> 0 —3 | | “ ss 
r Ss \ 
Evaluation of the determinants shows ¢(@) = 0 when a roll 
disturbance is applied. 
i 
b. @ due to yaw command (AB + By) (s) sae 
1 Je 
soyr Ss 
o{s)= —,—_ (2,16) 
i 
sal Ayr! ss 
Lim ¢(t)= Lim s-[,¢(s)] = = ee (2. 16a) 
t -co s~> 0 carmen] 
5° ss 


Again evaluation of the determinants shows ¢(®) = 0 when a yaw 


command is applied. 


ey ir due to a yaw command (AB+ By) (s)=—= 


IE 


ee 
Wb (s)= LS (271%) 
~ 
° ° P s° eee 
' _ 1; Se ee 
Lim  (t) = Lim s[ (s)] Al (2.17a) 


t-©o s>0 ss 


e 
Evaluation of the determinants shows that | goes to the constant 
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° -K, P 
[wy yee = K, cos @, + Ky (2.17b) 


and that J increases without limit as t~>™® 


e 
d. w due to roll disturbance 6p(s) = 


—— A 
is) (2.18) 
A 
. Rs°[A__] 
Lim  (t)= Lim s[w a ———— (2. 18a) 
t-o@ s > 0 [LA] .. 


Evaluation of the determinant shows that W (%) = 0, 
Further evaluation of the integral shows that | ©) = 0. 


This information on steady state values is of use in correlating 
the analogue computer results in@hapter 4 with the physical operation 


of the system. 
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CHAPTER 3 


COMPENSATION 


3,1 General Approach 


The general problem of preventing interaction between con- 
trolled quantities in multiple-loop systems has been discussed by 
Tsien, Pot Tires ‘ and others. Among these ena has developed 
a generalized theory of multiple-loop systems including techniques for 
the determination of compensation necessary to prevent interaction. 
Chu employs methods of manipulation with the system determinant 
with the aid of block diagrams, The methods of Tsien and Chu are 
quite parallel in theory, the choice between them being a subjective 
matter. Chu's method has been chosen for application to this 


problem since that technique is more familiar to the authors. 


A brief description of the basic elements of Chu's approach 
to mulitple-loop control systems is included in Appendix C and no 


further detailed explanation will be given here, 


The problem of analyzing the effects of interactions in this 
system is essentially one of maintaining contact with the physical 
hardware involved. Particularly, a block diagram possesses the 
general appeal of being a reasonably tangible representation of the 
relations among the physical variables of a system, provided that 
diagram does not endeavor to lump too many quantities and inter- 


relations into a single complicated block transfer function. 


A diagram which, it is believed, will provide a reasonably 
physical view of the problem will first be derived. It employs the 
most basic set of transfer functions obtainable from the linearized 


equations, 
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iP | at "i 
‘ie 


- 


Consider the set of linearized equations in the form ofa 


determinant: 

(I) (s-A) 

(II) -C 

(III) -H 

(IV) 0 

(V) -K,A 

where (I) 
(IT) 
(III) 
(IV) 
(V) 

and x 
Y 


From each 


cos oo -~sin a -B -~-M = 
S 0 -E -N = 
0 (s-F) -~L -G = 30 
x Z 
0 ou 0 l = 30 
S 
-Y 0 (s-K,B) -K|M =) 
(3.1) 


yaw force equation 


yaw moment equation 


roll moment equation 


roll control equation 


yaw control equation 


2 
OSE K,s + K,) 


(Ko + Ks) 


equation represented by a row of the determinant 


a single variable is selected to be described as a function of all the 


others. 


relation is defined: 





For manipulative ease, a transfer function describing each 


(3, 2) 


Equation (1) 


= F(s) (SieeZ ai) 


aoe 
WO 
i 
> 
~~” 
HI 


26 





sin a. 
(5) = MEZA) TA) = F.(s) (3, 2b) 


B . _2 - 
” > Toay = F.,(s) (3. Ze) 
<(s) 2 — = F,(s) (3, 2d) 


The theorem of superposition suggests representation of B(s) 


by the elemental block diagram. 


fn 
- 
oS 
Ss F + 
6 
S i, 
glen ot 
b B Y 5 b 
c lel ni o (SO) 
Equation (II) 
® (s) = Ls FAs) (3, 3a) 
| E _ 
wa ce F..(s) (3. 3) 
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0 
+ (s) = = = F(s) (3. 3c) 


@ 
+) ao (3, 3d) 





@ 
p B Y Ww 5 (3. 4) 


Equation (III) 


b Ho 

B ‘s) = (s-F) = F,(s) (3. 4a) 
£(s) = ra = F ,(s) (3. 4b) 
q G 

$ (s) = ess je F 9(s) (3. 4c) 
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$(s) 


6 
: Fo 
Fig. 3.3 








(3, 5) 

S(s)=—S = FCs) (3, 5a) 
S 

(3.6) 





Equation (V) 
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It will be noted that provision has been made for the insertion 
of a yaw command signal which subtracts from the sensed side 


acceleration (AB + By + M6) as indicated by Figure 1. 3 of Chapter lL. 


K,A 
ya eee). FZ) (3, 6a) 
BS) = SKB)? 
Y 
(s) a | (s) (3. 6b) 
$ (s-K,B) 13 
Y KUM 
88) © Te-KB) = Fil) (3. 6c) 
“Ky 
(AB+ By hm = (s-K By) = F(s) (3, 6d) 
B( 
: Fo 
4s) Ks) 
S S 
Fig ( } 


Fig. 3.5 


The interrelations existing among the system variables may 
then be described by a block diagram as shown in Figure 3.6, 
Provision of a node for introduction of a roll disturbance torque bp 


will be noted, 
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Fig, 3.6 System Block Diagram 





ol 





The method of Chu requires translation of the block diagram of 
a system such as depicted in Figure 3.6 into a "Standardized Block 
Diagram". In this "standardized" diagram certain manipulations to 
Figure 3,6 are made to satisfy a set of conventions which allow a 
symmetry to exist between the “standardized'' diagram and the system 
determinant and which allow the writing of the system determinant 
by inspection, These are described in Appendix C., When the 
necessary adjustments are made the "standardized" block diagram 


takes on the appearance of Figure 3.7. 


In Figure 3.7 it will be observed that the manner of redrawing 
Figure 3.6 to produce the "standardized block diagram" has resulted 
in the division of the coupled system into regions which can be referred 
to as the roll and yaw subsystems. It is readily apparent that Fo and 
Fy represent the coupling effects of yaw into roll while Fy, Fy, FY and 
Pi4 clearly show the coupling effects of roll into yaw. Inspection of 
the transfer functions represented by Fy, Fy, FL. Fy, Fy, Fig immedi- 
ately show that at oF = 0 the system separates into two subsystems, 


td 


By application of the rules (Appendix C) for writing the system 


determinant from the "standardized block diagram", that determinant 


becomes: 
Yaw Subsystem Roll to Yaw Couplings 
: =F ene | ee 
-Fe -F. | 0 -F, 
-F, -F, 1 | -F, -F, 
Pe gs (3.7) 
-Fy 0 -F, | ] “Fo 
0 0 0 | “Fu ] 
Yaw to Roll Couplings Roll Subsystem 
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Equating A to zero and evaluating one obtains the system 
characteristic equation. Immediate proof that above represents the 
original determinant of the system can be obtained by reading each 


row of A with the appropriate substitutions. 


As a consequence of the way selected for redrawing Figure 
3.6 to produce Figure 3.7, A divides into regions illustrative of the 
subsystems and the intercouplings as noted above. The close 
parallel between this determinant and Tsien's aie matrix are 
apparent. Also the clarity with which interactions are displayed in 
Figure 3.7 and in A would seem to some degree to refute one of the 


objections to the block diagram approach®" > 


It is true that considerable effort (although simple and direct) 
has been put into construction.of this form of A, but this has placed 
the designer in the position of having a graphic representation in 
Figure 3.7 of Cramer's Rule. Since each node output of Figure 3,7 
represents one of the variables of interest in the problem it is merely 
necessary to follow Cramer's Rule to assess the output at any node 


due to input at any node. 


3.2 Objectives and Procedures 


The objectives of complete compensation will be assumed to 
be: 
® 
IN isolation of ¢ from (AB + BY). 


2. isolation of UL from Onn 


The introduction of any compensation will change A, and in 
this will clearly affect all outputs. This change should not produce 
serious effect in the desired responses. (i.e. Should not produce 


an unacceptable array of closed loop poles). 


From the foregoing statements regarding the association 
between Cramer's Rule and Figure 3.7 and A it is clear that it is 


possible to solve for § (s) as a function of (AB + BY). as follows: 
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. F A 
any ore ey 


Where a has been formed from A by replacing the IV (or $) 
column of A by a column having zeros in all row elements except the 
first at which unity appears (the entrance function a and the yaw 
command signal transform having been factored out of the $ column of 


Ay). 


Likewise 


: A 
W anne 
oe aaa (3. 9) 


The desired response in yaw is given by 


A 
po a 
(AB+By)‘* > eee (3. 10) 


And the normal (existing, though not particularly desired) 


reaction of $ to Onn is 


° A 


_ re 
als) _— (3, 11) 





Provided the missile can be prevented from yawing due to roll 


disturbance, no side acceleration will occur, 


Then the desired goals of compensation are achieved if 


re ee = is = 0 
yr ry 
Ze aps Bp, is not altered by accomplishment 


of (1) except in the sense that A is changed, and 


that this change is not deleterious. 
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It is of interest to inspect both A and A 
yr iy: 


-—— +-,- oe 


13 
-Fe 1 -F, -F, 
renee 3) (UF) : b ~Fy (3,13) 
VAs 
“Fo 0 -F, 9 “Fig 
0 0 0 l 
| | p “Fig GO 8 -Fig 
“FP p 7h ote -P. 
ee ao b l Emr (3.14) 
-F 4 f -F, 1 -F 
+ -O-4 


As the dashed lines indicate, the expansion of ay and Da 
by cofactors pivots on the circled element and suppresses the 


indicated rows and columns. Thus 


l. acompensator placed at rJI, cIV will not appear in 


A 
yr 


2. acompensator placed atrV, clII will not appear in 
A 
ry 


Also of interest are A anaes fe 
cr yy 
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; 17 "i ae Va 
oy Fe f  -F, 

A= |-F, -F, 1 | -F, (3.15) 
-Fy 0 -F -F 


-Fe b -F, 0 “Fy 

Ae = uei, f 1 ote ey (3.16) 
-Fo | “EF, 1 a 
0 { 0 -F,, 1 


From these it is noted that a compensator placed at rl, cIlV does 


not appear in A a and a compensator placed atrV, cll does not appear 
on A ; 
yy 
Thus one concludes that compensation at rl, cIV can possibly 
be used to make A ee O and a compensator atrV, cll can possibly be 
used to make —T = 0 and that neither compensator will appear in 


A eueto| (A\ ee 
yy ite 


From the way in which the system determinant is constructed 
from the 'standardized block diagram"' it is apparent that the rI, clV 
element of A constitutes a feedback from node IV to node I in Figure 
3.7. Inspection of the diagram reveals that use of this path implies 
that $ is available and that modified é can be added algebraically to 

y. Likewise rV, cIl of 4 constitutes feed forward from node II to 
node V. i is physically tangible. The same question of ability to 


mix the modified | signal with 6 arises here, If the compensating 
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networks required are realizable either exactly or approximately this 


question of mixing can be studied further. 


Installation of compensation functions - V(s) at rV, cII and W(s) 
atrI, clV inA, where, conveniently, no element previously existed, 


results in the solutions, when ate and a ry are equated to zero: 


ae + (LA - BH)s + (HE - LC) cos a, 


(GE-LN)s + (-AGE-LMC - HNB +LAN+HME +GCB) oat 


V(s) = 


612s" + 143s - 146500 


Vis = “Bi. 25 +43200) (3.18) 
} - 94° l.2s + 42200) 


O 


a. ee + 47.8s - 1548 


vis} _ 19° I (100, 300)s + 37, 570] oe) 
O 


sind (s-A)[(CG-HMK, )s+K (-BCG-LAN-HME+LCM+AGE] 
Bs) = (s-K,B)[(GE-LN)s+(-AGE-LMC-NHB+LAN+HME+GCB)| 


s (s-F)[ Ns+(CM+K, ME-NA-BNK, )] 


1 


+ i SS SS PS SS SPSS PSS hs tg os pe 
(s-K,B)[(GE-LN) s+(-AGE-LMC-NHB+LAN+HME+GCB)] 


(3,20) 


Unfortunately neither of Tne eg ees is realizable with 


passive circuit elements exclusively. 


Upon examination of V(s) it is found that it can be divided into 








2 
Se 6. 71s 157s 1600 
a,=4 
-3 2 me = 
Pee ook Ss 478 x10 ~“s .. 15.5x10 
vis =P 7 Seon OC - “5 -ETGm | Teor 376 (3. 19a) 
fe) 
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An approximation using the last two terms of vis)] . is 
24 


employed in the analog computer simulation with the result that 
while addition of this approximation to V(s) does not create an 
unstable system, neither does it attenuate the roll disturbance due 
to a yaw command input. Thus the indication is that the first term 


is of dominant importance. 


Thus hope of a simple approximation to V(s) which would 
attenuate the unwanted roll response is lost. By arrangement of 
the analog simulator to provide an approximation to the first term of 


V(s) Fs above, about 30% decrease in the extremes of roll response 
2 4. 
to yaw command is achieved, indicating that at least partial decoupling 


can be achieved. Addition of the second term as above does not 

show any appreciable improvement over the first term alone. At 
this point available facilities on the computer are exhausted. The 
actual computer techniques employed are discussed in Chapter 4. 


Figure 5.7 shows the improvement achieved. 


Since, if isolation is to be achieved, simulation of the first 


term of V(s > iS of paramount importance, it would appear that a 
24 


; ¢ 
tandem pair of derivatives on the available \™ sensor output is 
required, In this connection it is worthy of note that the system can 
be driven unstable quite easily when more gain than specified is used 


in the double derivative term of V(s). 


The problem of introducing the modification to mix with 6 has 
not been examined here, Available literate. ss does not indicate 
any hope of introducing the compensating signal other than to 
differentiate it and mix it with é. rather than 6. That this procedure 


works has been verified on the analog computer. 


It will be noted in Figures 5.8 and 5.9 that the desired 
% 
responses to a yaw command, (A8+ By+M6) andw, have not been 


adversely affected by the approximate compensation employed. 


Relative to the possibility of using an expanded version of 


W(s) similar to that employed to approximate V(s), it will be noted 
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that W(s) is dependent on K), and thus the process of adjusting gains 
would involve adjustment of the compensator. This would result ina 
"tail-chase'' operation and is not felt to be practicable. Computer 
facilities available do not permit examination of the effects of 


approximate simulation of W(s). 


Since the desired approximate compensation for a 12° and 
2) 24° are widely different, it is apparent that if the compensation 
for a, = 24° is installed permanently, the compensation for a, = i 
will not be proper, even in an approximate sense, Also any 
compensation left in permanently will have the effect of creating 
coupling at oe 0. Then some form of logical operation is necessary. 
It is suggested that, since the a, = 0 and Oe ee responses are not 
greatly different (See Figures 5.1 and 5.2), the compensator for 


a 24° alone be installed and that this be switched in and out. 


Inasmuch as maneuvers which would be likely to develop 
extreme angles of attack would probably be presaged by some violent 
excursion of the seeker tracking line relative to the airframe 
longitudinal axis, it seems conceivable that this information could be 
used as a basis for the necessary compensation switching. This 


possibility has not been examined in detail. 


In conclusion, a method has been described for discussing 
compensation in a clear manner which involves explicit relation among 
the measured quantities available. A compensation which will 
theoretically completely isolate roll output from yaw command has 
been found and has been approximately verified on the analog 


Simulator. 
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CHAPTER 4 


COMPUTER SIMULATION 


4,1 Description of the Computer Setup 


The missile system is simulated on a Reeves Electronic 


Analog Computer (REAC). Outputs are recorded on a Sanborn 4 


channel recorder, 


The idealized system is instrumented, and can be 


represented by equations (4.1) through (4.5) below. These are 


rearrangements of equation (1.78) through (1. 82), and provide 5 


equations in 5 unknowns with each equation written to specify the 


highest derivative present of one of the variables, The functional 


layout of the equations is shown in Figure 4,1. * 


oe 


g 


re 


We 


e 


G(6+ 6.) + Fo + HB + Ly (4,1) 
-(Ky6 + Kd + Kg) (4, 2) 
(AB+ By+Mé) +sina $ - cosa b (4, 3) 
CBt+Eyt+N6 (4, 4) 


K,[(AB + By +M6)-(AB + By +MBé) |] 


3 (4, 5) 


+ Kb +K 


*Pigure D.l, Appendix D shows the actual scaled layout used. 
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The system of Figure 4,1 is the one on which all gain variation 
and adjustment tests are run. In addition, the compensator network of 


Figure 4,2, representing the approximate compensator, 


6. 71s" 1, 97s 


6 = = = 
A (s) = vis) 54° «84463 * +463 


(4,6) 








as developed in Section 3.2 is used in testing the effectiveness of 


decoupling compensation at a> 0. 


Some modifications of the layout of Appendix D are necessary 
to permit the insertion of the scaled compensator as a transfer function 
® 


from ¥ to 6. The modifications are minor, however, and the system 


is not shown in its modified form. The approximate derivative network 
4.1 


is from Johnson, Analog Computer Techniques. 









+ . 


Fig. 4.2 Analogue Computer Simulation of Decoupling 


Compensator 
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4,2 Gain Adjustment Techniques 
The following gains found by the methods of Section 2.2 to 


stabilize the system are used as a Starting set in adjustment to meet 


specifications. 
Kae (457) 
a 6.18 (4, 8) 
K, = oe (4, 9) 
K, - 27.9 (4,10) 
K, = Zo (4, 11) 
Ke = . 01985 (4.12) 


Roll adjustment is done by recording and observing ¢; yaw 


adjustments made on the recorded step response of (AB+ By + Mé). 


By observation of the recorded responses, cyclical gain 
adjustments in roll and yaw, and at the 3 pitch angles selected are 
made. It is found that K. and K, are to be set at or near their 


maxima. Next adjustments of K, and K, in roll, Ko and K) in yaw, 


are made to give specified stable sparen of the subsystems with 
Suitable damping of the coupled modes. Subsystem frequencies (and 
thus rise time) are dependent on K) and Ke; damping an Ko and K, ° 
These are the rough guides used in obtaining the final gains listed in 
Section 5.1. Section 5.2 shows system response to step 6,, and to 


ale 
step (AB + By+ M6) at the gains selected. * 


4.3 Compensator Simulation 


A scaled setup** of the compensator of Figure 4.2 is used in 
testing the possibility of eliminating or attenuating the effects of a yaw 


command on roll output. 





*Figure D,3, Appendix D shows correlation between analytically 
derived response of the uncoupled mode and computer response for 
a common set of gains. 


*x*kSee Figure D.2, Appendix D. 
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Responses of the system with and without compensator are 


compared in Figure 5.7 through 5. 9. 


45 





CHAPTER. 


RESULTS 
o.l1 Gains Selected 
Final gains selected are: 
K) = 32.0 (5. 1) 
Ky = 6. 40 (5. 2) 
ks oS 3 om (5. 3) 
SiS i203 (5. 4) 
K, = Zao (2770) 
Ke, = . 0199 (5. 6) 


0.2 System Step Responses 
Figures 5,1 through 5.9 are responses of ¢, (AB+By + M6), 


and wb to step inputs: 
bp = 1 degree (5. 7) 


(AB + By +M6) . = 1 deg/sec (5. 8) 
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0.3 System Response Parameters 


Reduction of response data from Figures 5.1 through 5.6 gives 


the uncompensated system response parameters show in Table 5.1. 


(AB +By + M68) ax 
deg 


sec 





Table 5.1 


Response Parameters of the Uncompensated System 





Table 5,2 is reduced data from Figures 5.7 through 5.9 and 


gives parameters of the compensated system at a, = an 


(AB.+By it MG) pee. 


deg/sec 





Table 5.2 


Response Parameters of the System as 
Compensated at 2 De 


*Roll and yaw subsystem time constants as defined in Section 1. 3. 
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CHAPTER 6 


CONCLUSIONS 


6.1 Achievement of Specifications 


The uncoupled response specifications of Section 1, 3 can be 
readily achieved. At ee 0, coupling exists which must be elimin- 
ated or damped to maintain stable, specified performance. Stable 
coupled response is obtained at the expense of more sluggish 


performance in the uncoupled mode. 


Satisfactory response in both coupled and uncoupled modes 
may be obtained, but requires computer simulation for practical 


investigation and setting of system; parameters. 


6.2 Elimination of Coupled Response by Compensation 


Decoupling by compensation is theoretically possible. However, 
the compensating devices cannot be achieved with purely passive 
networks. Also, some physical changes within the control systems 
are necessary to provide points at which compensating signals may be 
introduced and extracted. Third, since compensation is a function of 
angle of attack, there should be some manner of varying the compensator 


characteristics with angle of attack. 


Since some success has been achieved with approximations of 
theoretical compensators, a more intensive investigation of decoupling 
compensators might evolve a compensator which is both effective and 


physically practical. 


6.3 Effect of Non-Ideal Servos and Sensors 


Since the bandwidths of sensing devices and control surface 


O7 





servos commonly employed in missiles of this type are of the order of 
40 to 50 cps it is clear that the factors of (B. 25) are quite considerably 


removed from the poles (s-A) and (s-F) and from the origin. 


The results of factoring the ideal system characteristic equation 
indicate that its roots are not far removed from the poles thereof for the 
gains selected by computer study in Chapter 4. Then by inference, such 
additional roots as are introduced in the non-ideal system by the 
presence of control surface servos and sensing devices would be located 
relatively close to the poles introduced by those devices, and the 


residues ‘associated with these roots would be small. 


Therefore, mathematical representation by the idealized system 
is considered correct and no gross performance characteristics are 
thereby over looked. Additionally, as noted elsewhere, available 
analogue computer facilities do not permit simulation of the servos and 
sensing devices, Conversation with Professor W.L. Markey of M.I.T. 
relative to the relation of additional modes introduced by non-ideal 
instruments to the basic modes of such systems, (especially when these 
instruments have wide bandwidths compared to the natural modes of 
the system) would seem to comfirm that this simplifying action is in 


keeping with common engineering practice. 
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APPENDIX A 


AERODYNAMIC COEFFICIENTS 


Values of the aerodynamic coefficients defined by equations 
(1,53) through (1.62) and evaluated at the selected values of a, 
considered in this thesis are shown below in Table A.1. The values 
were obtained by Convair (Pomona) engineering personnel, Dimen- 


Sions are such as to make all angles considered to be measured in 


degrees. 

a = a. = ee a. = pas 
A - 373 - .282 - ,364 
B + .0705 + .0787 + .0787 
C +8. 42 -80.9 -119.5 
E -60.6 -69.0 - 69.0 
F ave 237 yee ed - 1,237 
G +1480 +1480 +1590 
H 0 -133, 5 +1269 
L 0 -133.5 -612 
M 0 -, 0149 -, 0526 
N 0 $13.3 +46. 8 

TABLE A. 1 


AERODYNAMIC COEFFICIENTS 
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APPENDIX B 


SYSTEM CHARACTERISTIC EQUATION 


B.1 Ideal System Characteristic Equation 


Expansion of the determinant (3.7) gives the ideal system 


characteristic equation in symbolic form 


A-1l-Figko Fs - FoF Fe - Fg Filo - Fl Fs - Fig kg 


+ Fol-Fg - Fig Fs Fo - Fo Fig + Fg Fig Fe! 


= es 


g F7Pgh3-F 


| eh ee eo 


Tee [FF g Fy 9 


nee -PO PAF 


107 F5F4 PAF e+Fig Fs Fs! 


+P l-Fig-Fyo Fa Fo- Fig Fg F3- FoF g- Fg FatFio Fig Fs! 


-F,, FF. - Feo eee, ee ee 


+ FP IFS Fig? Fyo Fo Fo- Fig Fg Fe t Fo Fs Fig- Fo Fat Fig Fic Fe! 
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Upon substituting the aerodynamic coefficients evaluated at 


a, = 24°, the C.E. becomes: 


ere (516 s” + (50 K, +65, 400 K,)s + 65,400 K 


1 2! 


- s“(stl. 237)[-s° +(-69K +.0787 K, - . 364) Pe - (103.5 K 


3 1 2) 


+ (Kg 3° + K, s + K,){-120, 000 s + 37, 200 K)] 


+s°(K, s° +K, s + K,){1590 s + (-92,8 K, + 511, 3)] 


] 


2 


+(K, 8° + K, s + K,)(K, s + K,)(81,200 s + 40,650) 


(B, 3) 


B.2 Non-Ideal system characteristic equation 


To obtain the characteristic equation of the coupled system 
when non-ideal sensing devices and non-ideal control surface servos 
are employed it is useful to note that these affect only the "roll control" 


and "yaw control" equations of (3.1). 


The "roll control" equation is modified to 
6(s) =-G.(K,8"+K-stK,)G, (45)¢(s) (B. 4) 
5° 6 5 4° 2 5° a 

The “yaw control'' equation is modified to 


sy(s)= G,[K,(AB +By+M6)G, + G(K,s + Ky)o] (B. 5) 


When these are case in the form employed in (3.1) the 


system determinant becomes 
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s-A cos @ - sing -B -M 
-C Ss 0 -E -N 
-H 0 s-F -L -G 
Go4G-xX 
0 0 4 0 1 
S 
(-K,G,G)A) (-G,G,Y) 0 (s-K)G,G)B) (-KG,G,M) 


(B. 6) 


where G, through G,. are as defined in Chapter 1Section1.2, This 


OD 
constitutes a modification of the transfer functions defined in Chapter 
3 based on these two equations, The following modified transfer 

functions are therefore defined to replace the corresponding transfer 


functions of Chapter 3. 


GaGEx 
s(e)=-—ye—= Fy, (B. 7) 
d S 
Y K)G,G A 
B oe (s-K,G,G,B) Pica (B. 8) 
Go4G,Y 
Y(s) = aes ee F a 
d (s-K,)G,G,B) l3a 
K,G ,G,M 
vies 0 olen Aaa ; 
ge ~ (s-K,G,G,B) ~ Fi4a (B, 10) 
1 4°71 
-K,G,G 
jae 
(AB+ By),'®) * (= G,GB) * Faa (B. 11) 
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With these modifications the same block diagrams, Figures 3.6 and 3.7, 
can be used and the expanded form of (3.7), see (B.1),be readily modified, 


For this purpose it is convenient to regroup the terms of (B.1) as follows: 


-1-F,, 1 -F.F, - F, - F 


3 6 He eee ake 


oF oF + FoFgks! 


+ Fiogl- FF, - Fg - FoFo! 


+ Pi aFisalFeF io -F5F4Fo Fy Fg F3-FoFy + FoF akg t FoF 5F3! 


+ If FLj-F F ,-F,F.F 


Mees 10> 6x4 iF 5F io! 


+ PuaFiaal “Fak + Fig Fg-Fl Fy FotF oF Fg! 


+ Fy oFiggl -FeFs - Fg - FFF, + FUFoFs! 


- PF. = FoF, = BOF. (B. 12) 


Upon substitution for all defined transfer functions except F 
through BY 


ll 


4 there results 
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Since the control surface servos and the sensors are presumed* 


to be second order devices, it is possible to assess the effedts of the 


modified transfer functions F,,_ through F as follows: 
lla l4a 
Define 
G = 2 — (B. 14) 
s +as +b 
Co Z (B. 15) 


soe 


Gs = (B. 16) 
s tes¢f 
G, = a= (B.17) 
Ss + §S 7A ee 
y n 
Ge = 5 (B. 18) 
S + Ms +n 
Then the modified transfer functions become 
Ee £ o - (B. 19) 


ee. ee 
lla (s +ces+d) (s +ms+n)s 


KA h b 


i : (s +esth) (s tas + ) (B. 20) 


$+ Ki? —5— eee ee 
(s A gst+h) (s +as+b) 


K,A hb 
ee (B, 21) 


[s (s° + gst hs? +as +b) - K,B hb] 


*Telephone conversation with Mr. R.L. Geisberg,CVAC, Pomona. 


66 





f h 


im ee 
se (meee =f) (aca sth 


s(s_ +gs+h)(s +as +b) - K,B hb 
(s° + gst A ee +as +b) 


: fh(s* +as +b) 


[ aie + gs t+ Ace +as +b) - K, Bhp] le “FSIS ap 18, 
(B. 22) 


KM + |. 


F Z (s tgsth) (s +as +b) 
l4a 


[s (eo + gs + pee +as+)- K, Bhp] 


(ee + gs + Ne +as +b) 


K, Mhb 
i. >. > (B, 24) 
Is(s +gsth)(s +as +h) - K, Bhb] 


From these it is clear that the common denominator of equation 


(B.13) will contain the product, 


s°(s°+cs+d)(s“+mgtn)[s(s°+gs+h)(s+as+b)-K, Bhb] (s“+es4f)(-Ae=F) 


(B. 24) 


Of these poles of A, ey Cen are as in the ideal system. The 


other poles are due to the factors of 


ee dy(a-mein)(s(s-+esth)(s-+as+b) - K,Bhb] eGaeee (B, 25) 


} 
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APPENDIX C 


CHU'S METHOD 


It is the purpose of this appendix to briefly describe the pro- 
cedures by which the ''system determinant'' employed in Chapters 2 
and 3 is derived from the “standard block diagram''. A complete 
development of these standardized manipulations as well as a con- 
siderable body of analysis and design theory for multiple loop 


systems can be found in the works of Chu: and Thelen oe Ay 


It is clear from the manipulation performed in the early part 
of Chapter 3 that the principle of block diagram representation of a 
set of linear integro-differential equations is simply one of expressing 
in LaPlace notation the relations in each equation as a summation 
process wherein several quantities are employed to form a single 
output quantity from anode. This process is repeated for each 
equation, and finally, interconnections between identical variables are 
made to complete the block diagram. This procedure is in common 


use. 


From the intermediate block diagram thus derived, the system 
characteristic equation can be gotten by any one of a number of 
methods including that due to Mason a8 as noted in Chapters 2 and 3, 
however, the method developed by Chu'©- eee additional facility 
in performing analysis for each of the several system variables and 


in design of compensation, 


68 





Figure C.1is the block diagram of an ''n node system" (n equations inn 


variables). 


The diagram permits of consideration of the introduction of 


any number of inputs (forcing functions) to any combination of nodes. Any 


block diagram derived by the common process stated above can be recast 


into the ''standardized'' form by simply listing the nodes from_left to 


right in any sequence and performing the simple diagram manipulations 


necessary to accommodate the following conventions: 


a. 


Each nodal point may receive any number of signals, but 


delivers only one output. 


Signals arriving at a node from a node to the left (feed 
forward) are arbitrarily considered to enter a node 
additively. (Sign adjustment in the feed forward transfer 


function may be necessary to accommodate this convention. ) 


Signals arriving at a node from a point to the right (feedback) 
are arbitrarily considered to enter a node subtractively. 
(Sign adjustment in the feedback transfer function may be 


necessary to accommodate this convention. ) 


Signals to be fed forward are derived at a connection 


immediately to the right of a node. 


Signals to be fed back are derived at a connection immedi- 
ately to the left of a node. (This and (d) frequently require 
moving a "pickoff point'' where a signal actually exists in 

the basic block diagram to the right or left past a transfer 


function block. ) 


Referring to Figure C.1 and considering the equilibrium equality 


at node a we have, for the sign conventions noted above: 


a=- (G_.G 


Ae (Gi -Gpa)h ise: (G_Gy,)n (C.1) 


+ G1 + Go, [i ae GiaW 
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This equation can be rewritten as: 


: m=V 
2S 
Clee G4Gnn)@ + (G, Gp.) ee F (Go Gna) ey (G_,,)m 
m= 
(C3) 
At node b we have the equality: 
b = (G_))a = (G, .Gpp) = (Go yGeap)e =e = (Go Gyp)? 
+G),1+ Gy I+..... cr (onr)) 
which can be rearranged to give: 
-(G_ Jat (it bcp?” +(G, Gaplet..-. +(G_ Gyp)n 
m=Y 
> (Gym (Ca) 
m= 


Extension of this process to several nodes using care to maintain the 
sequence established in (C.la) and (C. 2a) it can be seen that the 


sequence of nodal equations of equilibrium can be expressed in a 


determinant: 
PA Gye CBa Cod Ga 3 Gno'Na 
“Coys Dee Bb GoaSen «+--+ GnoGnp 
& = sore pone IHG aac no GONE (CLS) 
an or aoe HG, oeNn 


Inspection of (C.3) leads to the conclusion that, given the block diagram 
in standard form it is possible to write down the elements of the system 
determinant by inspection of the diagram, using only so much of the n 


node determinant pattern as applies to the system at hand. The main 


ie 





diagonal divides feedback and feed forward regions and terms of the 
main diagonal are the characteristic equations of the interior loops at 
the node indicated by the row and column of the main diagonal element. 
Expansion of the determinant in any given case produces the system 
characteristic equation without the often confusing mental bookkeeping 


necessary with some other methods, including that of Mason'@ a. 


By suitable choice of the method of selecting the sequence of 
nodes for drawing the "standard block diagram" it is possible (as in 
Chapter 3) to graphically illustrate in the determinant the couplings 
that exist within the original equations. Solution of the equations for 
desired input - output relations follows Cramer's Rule, the column 
into which the forcing functions are substituted, representing the node 
in the standard block diagram at which the desired variable exists. 


The system being linear, superposition is applicable, 
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APPENDIX D 


SCALED COMPUTER SETUP AND CORRELATION RESULTS 


D.l. Scaled Layout 


The REAC layout of the simulated system is shown in Figure 
D.1. Note that a time scale factor of 10:1 is used (i.e. , response is 
slowed so that the computer takes 10 seconds to respond as the real 


system does in 1 second). 


D.2. Scaled Compensator 


The compensator V(s) of Chapter 3 is approximated as shown 
in Figure D.2. It is limited in both the inability to produce a perfect 
derivative and by some approximation of time constants to conform to 


the RC values available, 
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Fig. D.2 Scaled Compensator Layout 


73 





(s) 


| seps008d era FePTOT 


sca. = 


I 
O1 
O1 


(st pen - 


*¢ 


Jap10s3ay 
a 
| 





wa}Shg pazijeap] jo NOAeTY DYaY PeTeOS 1a std 


(19 4QOT 


1 . - 
Yon+ Lat+dy) oe 


2 ©) 
7 
© a << ; Ol or (91) ees 
pa: 13 
; @ a Co cmars V)0T- 


(orerLat ey) yS- 





c ar 2a 
| or ee ql - A @ Ag- 
“T ell cD 


T OTs 


PINOT (0% 9OI- , 
cVvoOl (81) 
HO! G S 


0 
1 - 


D SOdn - w i 
- o oS 
a . 


(SIN +*98 + E9V)oI- 


Vo UIS 


goI- gol 


[pe= 2] © 


Th owt+*at gy)- 


Go. 
nN 
— 
" 
° 
o Ss 


1o® = 
@) Ar 08 - S eS 
y gee (01) Hoe © Ol; 
o y I o§- ! 02 Lg 50%. (e) A00T- 
£5 
Y gf B0c. = 90I- 
oT ~ 


*-o 
ae 
—|S 
1 
*- 
aN 


74 





D.3 Response Correlation 


Figures D. 3 and D.4 are correlation check runs between 
computer and analytic responses of the uncoupled roll and yaw sub- 


systems. 


Gains used are not the final selected values of equations (5. 1) 
through (5.6) but are listed below: 


kn eee ae (D. 1) 
Say 6.18 (1922) 
Lea) = noon (D. 3) 
K, = 27.9 (D. 4) 
nae 2 Oe (D. 5) 
Lor , 01985 (D. 6) 


When these values, analytic responses below are obtained: 





mimeeeiptede--- O's 4g56 29S (568 2a ee am 
(D. 7) 
and 
ieeriesos 79ee '-°' + 907e (°° ain(1G, Gareinomeney 
(D. 8) 
where 
7 _ , deg 
P= |(AB+ By ), = 1 So8 (D. 9) 
and 
R = |S = 1 deg (D. 10) 


In Figures D. 3 and D.4, the continuous lines are computer 


response, circled points are check points from equations(D.7 and (D. 8). 
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